Fractal Conductance Fluctuations in a Soft Wall Stadium and a Sinai Billiard 
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Conductance fluctuations have been studied in a soft wall stadium and a Sinai billiard defined 
by electrostatic gates on a high mobility semiconductor heterojunction. These reproducible mag- 
netoconductance fluctuations are found to be fractal confirming recent theoretical predictions of 
quantum signatures in classically mixed (regular and chaotic) systems. The fractal character of 
the fluctuations provides direct evidence for a hierarchical phase space structure at the boundary 
between regular and chaotic motion. 
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Ballistic geometries defined in two dimensional elec- 
tron gases (2DEGs) provide excellent systems for study- 
ing quantum chaotic phenomena. In magnetic fields these 
systems yield reproducible conductance fluctuations and 
weak localization effects analagous to the universal con- 
ductance fluctuations and weak localization effects stud- 
ied in disordered conductors. In the ballistic case the 
phase coherent phenomena reflect the device geometry; 
not the random positions of impurities. By studying 
these effects in both chaotic and non-chaotic geometries 
a valuable insight has been obtained on how quantum 
information can be retrieved from classical chaotic dy- 
namics. To date, in virtually all of the systems studied 
both experimentally and theoretically, a purely chaotic 
system has been assumed in which classical trajectories 
probe the phase space in a purely ergodic fashion 
before they exit. Escape from the ballistic cavity in such 
systems usually occurs exponentially fast. 

Real billiards, however, are typically not fully chaotic. 
They have a mixed phase space, i.e., they contain both 
chaotic and regular regions. Naively, one might expect 
that they are simply combinations of independent fully 
chaotic and regular regions. But this is not the case. 
Mixed systems possess an important property which has 
drastic consequences for conductance fluctuations: the 
chaotic part of phase space obeys a power law escape 
probability Q, in contrast to the much faster exponen- 
tial decay of fully chaotic systems. The power law origi- 
nates from chaotic trajectories which are 'trapped' close 
to an infinite hierarchy of regular regions at the bound- 
ary between the regular and chaotic motion (Fig. 3). This 
trapping is believed to be a consequence of Cantori, which 
act as partial barriers for transport and which surround 
the regular regions at all levels of the hierarchy HI). 
The phase of these long trapped trajectories is extremely 
sensitive to any externally changed parameter. In a 
magnetic field, for example, they acquire a phase factor 
exp(2niAB /cf>o), which depends on the number AB/cpQ 
of magnetic flux quanta 4>q = h/e enclosed by the trajec- 
tory, (A is the accumulated enclosed area). The larger 



the value of A, the smaller the change in magnetic field re- 
quired to modify the phase. Therefore, in these systems, 
one expects to observe conductance fluctuations with 
much finer field scales |f7U|. Semiclassical calculations 
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FIG. 1. Resistance vs. magnetic field (a) for the open stadium 
at a gate voltage of -1.9V after illumination at 50mK and (b) for the 
Sinai Billiard at 50mK. Fluctuations on both large and small scales 
can be seen for both devices. The two insets are SEM micrographs 
of devices similar to the ones used. In the actual Sinai billiard used 
for these measurements an insulator bridge connection (not shown) 
was used to make contact to the center gate of the Sinai device. 
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similar to those in Refs. have demonstrated that 

a power-law distribution of enclosed areas larger than A, 

P(A) ~ A-\ (1) 

leads to fractal conductance fluctuations [|| . Under these 
conditions the line conductance versus magnetic field has 
the same properties as fractional Brownian motion 
which is self-affine in a statistical sense and which is de- 
scribed by a fractal dimension, D, given by 

D = 2-7/2. (2) 

Thus the exponent 7, which describes a property of the 
classical phase space, is related to the fractal dimension 
D of a quantum coherent measurement. Similarly, the 
variance of conductance increments scale with small mag- 
netic field increments as < (AG) 2 > ~ (AB) 1 . These 
results hold || for 7 < 2. Since 7 > 1 must hold from 
more general arguments |l0| this restricts the fractal di- 
mension to lie between 1 (7 = 2) and 1.5 (7 = 1). While 
the value of the classical exponent 7 is non-universal and 
is found numerically to be sensitive to details of the ge- 
ometry and the confining potential S, the occurrence of 
power-law distributions is universal || . It is not feasable 
to test the existence of fractal conductance fluctuations 
by numerical quantum calculations, since a large number 
of modes are needed to adequately probe the hierarchical 
classical phase space. Indeed, because the semiclassical 
calculations of Ref. Q are based on many assumptions 
and approximations, one may even speculate if fractal 
fluctuations exist at all. 

The first experimental evidence for fractal behaviour 
was recently found in gold nanowires p]. The variance 
analysis of the fluctuations showed a power law behav- 
ior over one order of magnitude in magnetic field. While 
suggestive, this experiment suffered from two limitations, 
i) it is well known that almost any smooth function can 
be satisfactorily fitted by a power law for just one or- 
der of magnitude, and ii) for the rectangular geometry of 
these gold nanowires one does not expect a mixed phase 
space, unless fortuitous fluctuations in wire width occa- 
sionally create such regions |M. Herein we report the 
first observation of fractal conductance fluctuations over 
two orders of magnitude in magnetic field and in genuine 
mixed phase space geometries. The two classic models 
of fully chaotic systems, namely the stadium (two half 
circles connected by straight lines) and the Sinai billiard 
(a square with a circular disk at its center) are defined 
by electrostatic gates on a high mobility semiconductor 
heterojunction. These systems are no longer fully chaotic 
due to soft-wall potentials, but have a mixed phase space. 
This can be easily verified numerically (Fig. 3). We find 
fractal conductance fluctuations as a function of mag- 
netic field for large ranges of applied gate voltage (corre- 
sponding effectively to different billiard sizes and poten- 
tial forms). The existence of this new quantum signature 



of classically mixed systems is therefore conclusively con- 
firmed. In addition, by choosing the stadium geometry, 
we have shown that the famous experiment by Marcus 
et al. on the stadium ||, which was performed to study 
the properties of fully chaotic systems, gives qualitatively 
different results when carried out on todays high mobility 
samples. 

The two devices are shown as insets in Fig. 1. The sta- 
dium had a lithographic radius of 1.1/im. It was defined 
using metallic gates on a high mobility AlGaAs/GaAs 
wafer (with a mobility of 2.2 x 10 6 and 4.3 x I0 6 cm 2 /Vs 
and density of 1.7 x 10 11 and 3.3 x 10 11 cm before and 
after illumination with a red LED). The 2DEG was 95nm 
below the surface. The device leads were made unusually 
wide (0.7 microns). This had the effect of allowing most 
trajectories to rapidly exit the stadium with the excep- 
tion of those trajectories which were trapped near the 
hierarchical phase space structure at the boundary be- 
tween regular and chaotic motion. This feature and the 
high mobility wafer used made this an optimal device for 
the observation of fractal conductance fluctuations. The 
very high mobility wafer was necessary to achieve the 
required long phase coherence length, l v . The sub mT 
features in the conductance fluctuations and the narrow 
weak localization peaks (< than 150/iT) confirmed that 
Itp was indeed many times the perimeter of the stadium 
and Sinai billiard. Lithographic and fabrication details 
of the Sinai billard have been published elsewhere [jl3| . 
The measurements were made on a dilution refrigerator 
using standard low power AC techniques at 50mK. 

Measurements at higher temperatures (~ AK ) revealed 
small features in the magnetoresistance related to the 
classical focusing of trajectories. As the temperature was 
lowered both conductance fluctuations and a weak local- 
ization peak (around B=0) developed. In this paper we 
analyze the fluctuations which occur at magnetic fields 
for which the cyclotron diameter is larger than the sta- 
dium diameter. Features related to noise can be elimi- 
nated by comparing the ±B traces. In Ref. [Jl3| it was 
observed for a soft wall Sinai billiard that by rescaling 
magnetic field and resistance one finds quite similar se- 
quences of maxima and minima. This non- statistical self- 
affinity remains unexplained so far and does not occur for 
the stadium, suggesting that this feature is characteris- 
tic of the Sinai billiard geometry. We observe in contrast 
conductance fluctuations on very different scales resem- 
bling the statistical self-affinity of fractional Brownian 
motion for both geometries and for large ranges of gate 
voltage. 

Experimental traces for both geometries are shown in 
Fig. 1. The results of a fractal analysis on these curves, 
shown in Fig. 2, gives fractal dimensions of D = 1.25 
and D = 1.30. The fluctuations are found to obey 
power-law scaling for over two orders of magnitude in 
magnetic field. In all observable fractals this scaling be- 
haviour is bounded by upper and lower cutoffs. In our 
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case the upper cutoff, in magnetic field, is determined 
by the smallest area for which the power law distribu- 
tion holds (which is close to the range of magnetic fields 
experimentally studied). The lower cutoff is determined 
by the minimum of two time scales: i) the finite phase 
coherence time and ii) the Heisenberg time tn — K/AE, 
when the average level spacing AE of the closed device 
is resolved and the semiclassical approximation becomes 
unreliable. The power law area distribution for trajecto- 
ries staying in the device longer than any of these times 
will not lead to fractal fluctuations. The fractal dimen- 
sion observed between these cutoffs is determined by us- 
ing a refined version of the box-counting algorithm. In 
the standard box-counting algorithm one puts a grid of 
square boxes of size LxLon the data (conductance vs. 
magnetic field) and counts the number N(L) of boxes 
through which the curve passes. Its dependence on box 
size L, N(L) ~ L~ D , defines the fractal dimension D. In 
this standard analysis the relative scale of conductance 
and magnetic field is arbitrary so that for a finite data set 
the resulting fractal dimension depends on the aspect ra- 
tio of the plot. To overcome this limitation we applied a 
modified version of the algorithm: This divides the mag- 
netic field range in length intervals AB and determines 
N(AB) as the difference of maximum and minimum con- 
ductance in each interval, summed over all intervals, and 
divided by AB. This corresponds in the standard algo- 
rithm to taking rectangular boxes with infinitely small 
size in the conductance direction. The inset of Fig. 2 
shows the analysis of the variance, which gives 7 = 1.38 
and 7 = 1.30, for stadium and Sinai billiard respectively, 
confirming the above fractal analysis JT1]. 
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FIG. 2. Fractal analysis for the stadium (diamonds) and 
Sinai device (triangles) data shown in Fig. 1 using a modified 
box-counting algorithm. The number of boxes N follows power 
laws (AB)~ D for two ordes of magnitude in magnetic field scale, 
giving fractal dimensions D = 1.25 and D = 1.30, resp. The inset 
shows the variance analysis giving 7 = 1.38 and 7 = 1.30, resp., in 
reasonable agreement with Eq. (2). 
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FIG. 3. (a) Poincare surface of section of the stadium of Fig. 1 
(with OA/im depletion length and a OAfim wide parabolic wall). It 
shows v x vs. x at every trajectory intersection with the horizontal 
symmetry line whenever v y > holds, for one chaotic and 6 regu- 
lar trajectories, (b) Enlargement of the small island in (a) showing 
higher order islands, (c) Periodic trajectories corresponding to the 
center of the large island in (a) (dashed line) and the center of the 
island in (b) (solid line), (d) Quasiperiodic trajectory correspond- 
ing to the island in (b). (e) Chaotic trajectory being trapped in 
the stadium close to the regular trajectory of (d). (f) Integrated 
area distribution displaying a power law for large enclosed areas. 

We now show that the stadium, a classic model for fully 
chaotic systems, has a mixed classical phase space when 
achieved electrostatically due to the soft-wall potential 
experienced by the electrons. We model this potential 
by a depletion region of lOOnm at the edge of the device 
followed by a 400nm parabolic region and a flat potential 
beyond that The mixed phase space is revealed 

in a Poincare section analysis (Fig. 3). Its hierarchical 
structure gives rise to long trapped trajectories as well 
as a power-law area distribution. The effective Planck's 
constant % e g = %/ J pdq — (L\/2irn)~ 1 for a trajectory 
along the circumference of the stadium is fi e s = 8.6 x 
1CP 4 . All of this confirms the applicability of the results 
of Ref. H to the present experiment. Similar conclusions 
apply to the soft-wall Sinai billiard. 
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Direct Simulation is the only technique available to 
predict the exact value of the classical power-law expo- 
nent and the corresponding fractal dimension for a par- 
ticular device. The values of these exponents depend 
critically on the exact form of the potential || . But even 
if the electrostatic potential due to the gates is known, 
it would be modified in a real device by the presence of 
disorder. In high-mobility samples, however, the disor- 
der potential is weak and smooth and so does not change 
the qualitative character of the classical phase space, i.e., 
it remains a mixed phase space. The exponent 7 of the 
power law area distribution of trapped chaotic trajecto- 
ries on the other hand, might well be changed by the 
specific disorder configuration inside the cavity. Thus 
one cannot hope for a quantitative comparison of the 
exponent 7 deduced from classical simulations without 
disorder with the fractal dimension D of the conduc- 
tance fluctuations of a real device containing disorder. 
Disorder outside the cavity may in principle also change 
the conductance fluctuations but it cannot influence the 
trapping of chaotic trajectories inside the cavity and thus 
it will have no significance for the observed fractal di- 
mension. As a function of the applied gate voltage, the 
exponents of the classical power law distribution may 
fluctuate for large times |ll| but for times smaller than 
the phase coherence time the power law exponents have 
been found to be much more stable This feature is 
also reproduced in our experiment. Figure 4 reveals that 
there exists a monotonic relationship between the fractal 
dimension and the gate voltage applied to the stadium 
device. The fractal dimension tends towards 1 as the 
confinement is reduced. 



1.6 

2-9.8 

1.5- is. 9 - 4 

o 9.0 



-1.9V -1.5V 
16 



1.4- 
Q 1.3- 
1.2- 
1.1 - 
1.0 



-50 50 -50 50 -50 50 

B (mT) 



-2.0 



-1. 



-1.0 



-1.6 -1.4 -1.2 

gate voltage (V) 

FIG. 4. Fractal dimension of conductance fluctuations vs. gate 
voltage for the stadium before illumination. Error bars of the frac- 
tal analysis are typically ±0.05. The insets show data for three 
gate voltages. 

We have reported the first observation of fractal con- 
ductance fluctuations in semiconductor billiard devices 
fl7| ]. The observation confirms a recent theoretical pre- 
diction that conductance fluctuations in mixed phase 
space systems are statistical self-affine and can be de- 
scribed by a fractal dimension. The origin of the be- 



haviour is a power law distribution of areas enclosed by 
chaotic trajectories, which results from the hierarchical 
structure of phase space at the boundary of regular and 
chaotic motion. We have shown that in real devices the 
classification into 'chaotic' and 'regular' geometries is in- 
complete at best. Our findings confirm the important 
role that soft-wall potentials play in nanostructures. We 
have also observed a dependence of the fractal dimension 
of conductance fluctuations on the device parameters. 
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